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ABSTRACT 

The aim of this paper is to clarify the notion and cause of overmerging in N-body 
simulations, and to present analytical estimates for its timescale. Overmerging is the 
disruption of subhaloes within embedding haloes due to numerical problems connected 
with the discreteness of N-body dynamics. It is shown that the process responsible for 
overmerging is particle-subhalo two-body heating. Various solutions to the overmerg- 
ing problem are discussed. 

Key words: galaxies: evolution - dark matter - large-scale structure of Universe 



> 

(N 
O 
(N 
O 
O 

o 

S3 

Oh' 

6 

CO 



X 



1 INTRODUCTION 

In the study of the formation and evolution of large-scale 
structure in the universe, clusters of galaxies, galaxies, and 
many other gravitational systems, haloes play an important 
role. A halo is defined as a collapsed and virialized den- 
sity maximum. A halo can contain several smaller haloes, 
denoted as subhaloes. We denote a halo that contains sub- 
haloes as an embedding halo. The existence of subhaloes 
is an important issue in cosmology, especially for galaxies, 
and groups and clusters of galaxies. For example, hierar- 
chical structure formation scenarios for CDM-like spectra 
predict many more dwarf galaxies (or satellites) than ob- 
served (Klypin et al. 1999b, Moore et al. 1999). Also, a spiral 
galaxy cannot retain its disk if there is an high abundance 
of subhaloes within its halo (Moore et al. 1999). Thus, the 
question is whether the initial density fluctuation spetrum is 
such that not many dwarf galaxies form in the first place, or 
that they are easily destroyed within our Galaxy, and hard 
to find outside it. On larger scale, clusters of galaxies do 
contain an abundance of subhaloes, i.e. its member galax- 
ies, which are often distorted and stripped, but probably not 
destroyed. 

N-body simulations are routinely used to model the for- 
mation, evolution, and clustering of galaxy and galaxy clus- 
ter haloes. However, the N-body simulation method does 
have its limitiations, and care should be taken with the in- 
terpretation of the simulation results. One such limitation 
arises from the use of particles to represent the mass dis- 
tribution whose evolution one tries to simulate. If the phys- 
ical mass distribution is effectively collisionless, as is often 
the case, the use of particles gives rise to artificial colli- 
sional effects within the numerical mass distribution, espe- 
cially though two-body interaction. Two-body encounters 
between simulation particles, either close or distant, deflect 
their orbits significantly, while they should behave like test 
particles, and respond only to the mean potential. This is es- 
pecially a problem for subhaloes, which are easily destroyed 



in an N-body simulation with insufficient resolution (White 
et al. 1987; Carlberg 1994; van Kampen 1995). We use the 
term overmerging for the numerical processes that artifi- 
cially merge haloes and subhaloes in an N-body model, usu- 
ally by disrupting the subhalo. Thus, the term merging only 
denotes merging due to physical processes. 

However, there seems to be some confusion in the lit- 
erature over the nature, cause, and importance of the over- 
merging problem. This paper attempts to clarify the differ- 
ence between the three most important two-body processes 
operating on subhaloes, and provide estimates for their as- 
sociated timescales. The first process is two-body evapora- 
tion, which is an internal process operation within any halo 
or subhalo. It is due to two-body interactions between the 
particles within the halo or subhalo. The second process is 
particle-subhalo two-body heating, which is the heating of 
'cold' subhaloes by particles from the 'hot' embedding halo 
through two-body interactions. The third process is particle- 
subhalo tidal heating, where the subhalo is considered colli- 
sionless, and increases its kinetic energy through tidal inter- 
actions with particles from the embedding halo. All three 
process are illustrated in Fig. 1. 

Besides the two-body processes, the use of softened par- 
ticles, in order to minimize two-body effects, can cause over- 
merging as well by artificially enhancing physical processes 
like merging and disruption by tidal forces (van Kampen 
1995; Moore et al. 1996). Groups of softened particles are 
not as compact as real haloes, and their artificially larger 
sizes make N-body groups more prone to tidal disruption. 

Carlberg (1994) proposed particle-subhalo two-body 
heating as the main cause for overmerging. He gave a 
timescale for this process, but no derivation. This was pro- 
vided by van Kampen (1995). It has been shown before 
(Carlberg 1994; van Kampen 1995) that the two-body heat- 
ing time-scale for small subhaloes orbiting an embedding 
halo is short enough to result in their complete destruction 
and dispersion. Subsequent authors, including Moore et al. 
(1996) and Klypin et al. (1999a), referenced Carlberg (1994) 
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Figure 1. Graphical illustration of the main numerical disrup- 
tion processes that cause overmerging. Open circles represent the 
'cold' particles of an isolated halo or a subhalo within an embed- 
ded halo, whose 'hot' particles are indicated by filled circles. 

as saying that 'particle-halo heating' is at the root of the 
problem. Moore et al. (1996) then claim that the process is 
not important for the resolution of the simulation performed 
by Carlberg (1994) because the timescale is too long. How- 
ever, the process Moore et al. (1996) actually describes and 
derives a timescale for is a different one, driven by tidal 
encounters between simulation particles and perfectly colli- 
sionless subhaloes, while Carlberg (1994) and van Kampen 
(1995) clearly had a collisonal process in mind, driven by 
two-body encounters between individual subhalo particles 
and particles from the embedding halo. This paper shows 
that the latter process has a much shorter timescale, and 
therefore, along with excessive softening, is the main cause 
for overmerging. 



2 TWO-BODY EFFECTS IN N-BODY 
SIMULATIONS 

Two-body effects become dominant for systems modelled by 
small numbers of particles. This is usually quantified by the 
two-body relaxation timescale, which is defined as the time 
it takes, on average, for a particle to change its velocity by 
of order itself. After this time a system is denoted as relaxed. 
The relaxation timescale is defined in terms of the half-mass 
radius r, the typical velocity v (usually taken to be equal 
to the velocity dispersion), and the average change Av per 
crossing time r/v: 



tr 



(Av) 2 V 



(1) 



For an isolated virialized system of N point particles, it is 
easy to show that this is of order O.lJV/lniV crossing times 
(eg. Binney & Tremaine 1987). If the time interval one tries 
to cover for a particular problem is larger than the two- 
body relaxtion timescale, the problem becomes artificially 
collisional. 

Although many systems are likely to endure physical 
mechanisms like violent relaxation and phase mixing dur- 
ing some stage of their evolution, an obvious worry is that 



such physical mechanisms might not actually be important 
in a given situation, so that two-body interactions are very 
much unwanted as they might mimic the effects of physi- 
cal mechanisms. A different problem is that a system might 
completely evaporate through two-body interactions. 

The problem can be alleviated somewhat by soften- 
ing the particles, which reduces the two-body relaxation 
timescale to OTA/ In A crossing times (see van Kampen 
(1995) and references theirein), where A = Min(i?/4e, N), 
with R the effective size of the system, which we take to be 
twice the half-mass radius r, and e the softening length of the 
N-body particles. However, e is necessarily a function of both 
N and r, as softened particles should not overlap too much. 
Too large a choice for e will prevent particles from clus- 
tering properly and produce haloes which are too extended 
and too cold (that is, the velocity dispersion is too small). 
Given that N/2 particles reside within r by definition, the 
mean particle number density within r is 3N/(8irr 3 ). The 
maximum mean particle density desirable is set by the min- 
imum mean nearest neighbour distance for the particles: 
n. m ax = 3/ (4:TYr'l n ) . Most often used is Plummer softening, 
which just means that particles have a Plummer density 
profile, p(r) ~ (r 2 + e 2 ) 5 / 2 . The effective force resolution, 
defined as the separation between two particles for which 
the radial component of the softened force between them is 
half its Newtonian value, is « 2.6e for Plummer softening 
(Gelb & Bertschinger 1994), so we want r nn Jj 2.6e, which 
gives n max « 0.014/e . Thus, we find a maximum realistic 
softening length 

(2) 
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For this e the relaxation time becomes 0.3A/lnA crossing 
times, i.e. three times larger than for the point particles case. 

Even though softening alleviates the problem of two- 
body effects somewhat, softened particle groups are more 
extended and less strongly bound (van Kampen 1995). This 
makes them more vulnerable to two-body disruption pro- 
cesses, which are more efficient for larger subhaloes, as 
shown below. Furthermore, the timescales for physical dis- 
ruption processes are effected. Subhalo-subhalo tidal heating 
has a timescale inversely proportional to the subhalo size 
(van Kampen 2000), and is therefore slower, although the 
lower binding energy might compensate for this. Tidal strip- 
ping and disruption will be artificially enhanced, however, 
because of the larger subhalo size and the weaker binding of 
the particles inside the group (van Kampen 1995). Because 
the enhanced tidal disruption due to softening has the same 
net effect as two-body disruption, which is also enhanced 
due to the larger subhalo size, the two disruption processes 
accelerate each other. In the next section we derive two-body 
disruption timescales without taking into account tidal dis- 
ruption, and then treat these timescales as upper limits. 



3 DISRUPTION TIMESCALES 
3.1 Two-body evaporation 

This process is internal to haloes and subhaloes in other 
words, it is a self-disruption process. Two-body interactions 
between particles within the same (sub) halo change their 
orbits and velocities, thus every once in a while the velocity 
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will be larger than the escape velocity, and a particle will 
'leak' out of the (sub)halo. The timescale for this process is 
about a hundred times the relaxation timescale (e.g. Binney 
& Tremaine 1987), 

N r 

t dis = 30—-, (3) 
In N v 

so for small N this becomes important. As an example, for 
galaxy haloes, evaporation becomes an issue for N < 10, as 
the crossing time for most galaxy haloes is larger than 0.2 
Gyr, independent of their mass. 

Moore et al. (1996) tested whether this process gets en- 
hanced for subhaloes within embedding haloes due to the in- 
fluence of the mean tidal field of the embedding halo. They 
simulated a collisional group of particles within a smooth, 
and therefore collisionless, isothermal system. They found 
that the evaporation rate was similar to that for an isolated 
group, and concluded that "relaxation effects are not impor- 
tant at driving mass loss from haloes within current simu- 
lations". Their conclusion is incorrect, however, as they did 
not consider the particle-subhalo two-body heating process, 
which we discuss next. 

3.2 Particle-subhalo two-body heating 

Particles within a subhalo do not just interact amongst 
themselves (driving the evaporation process describe above), 
but also with the particles of the embedding halo. As the 
latter are usually hotter than those of the subhalo, veloc- 
ity changes to the subhalo particles will always be positive. 
The process very much resembles the kinetic heating of a 
cold system that is introduced into a hot bath: an embed- 
ding halo 'boils' the subhalo into dissolution. A derivation 
for the disruption timescale of this process is given by van 
Kampen (1995, his eq. 15, which is erroneous by a factor of 
two): 
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Here TV denotes the number of particles, and the subscripts 
h and s denote embedding halo and subhalo respectively. A 
similar expression was given earlier by Carlberg (1994, his 
eq. 13 with his indices c and g swapped, no derivation given): 
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We can rewrite eq. (4), using eq. (2) and the virial theorem 
for both halo and subhalo, as 



t d 
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(6) 



4 In iV h r B w h 8 In JV h 

This timescale is shorter than that for two-body evapora- 
tion, by a factor of (van Kampen 1995) 

Vh In N h 



100 



v s rl In N s ' 



(7) 



where the virial theorem, v 2 ~ N/r, is used for both systems. 

Because r h is at least several times r B , the disruption 
time (6) is at least « 7V s /ln(7Vh) embedding halo crossing 
times, which covers the range 0.05 — 0.15iV s crossing times 
for iV h « 10 3 - 10 9 . Thus, it is a much faster process than 
two-body evaporation. 



3.3 Particle-subhalo tidal heating 

A different cause for overmerging was proposed by Moore 
et al. (1996): the tidal heating of subhaloes by particles of 
their embedding haloes. Subhaloes are taken to be collision- 
less, and get disrupted through an increase of their internal 
kinetic energy by tidal distortion from passing N-body par- 
ticles, which are artificially large as compared to the true 
dark matter halo particles. 

The time-scale for this process as given by Moore et al. 
(1996; their eq. (3), which is eq. (7-67) of Binney & Tremaine 
1987 with the assumption that the r.m.s. radius is equal to 
the half-mass radius) reads 



t dis « 0.03 



Vh m s r p 
Gn p m 2 r'i 



(8) 



where the subscript p stands for perturber. The perturber is 
an N-body particle of the embedding halo with mass m p 
and size r p , at a distance q from the centre of the em- 
bedding halo. Note that the impulse approximation implies 
v 2 /Av 2 = E/AE. Moore et al. (1996) then assume the em- 
bedding halo to be isothermal, so that n p « vl/2-KGm v q 2 , 
set the half-mass radius of the subhalo equal to the tidal 
radius, qv s /(3vh), and assume the subhalo to be virialized. 
This gives 



(9) 



Relation (8) was originally derived by Spitzer (1958) for 
the disruption by giant molecular clouds of open star clus- 
ters. An important assumption in its derivation is the tidal 
approximation, which is only valid for impact parameters 
b > 6min- Aguilar & White (1985) found that 6 m i n should 
be at least five times the size of both the perturber and 
the perturbed system. Binney & Tremaine (1987) use the 
tidal approximation down to brain — ^cluster < ^cloud , and 
introduce a correction factor g — 3 to take into account the 
encounters for which the tidal approximation fails. 

Moore et al. (1996) take this result and apply it to 
tidal interactions between the N-body particles of an em- 
bedding halo and its subhaloes theirin. Thus, they set r p to 
the gravitational softening length e. However, as the size of 
the perturbers is now smaller than the size of the perturbed 
subhaloes, the tidal approximation, even with the correc- 
tion factor g included, is only valid for b > r s . Therefore, 
setting r p = e in eq. (9) is incorrect; instead, one should 
set r p = r s . This means that the timescale becomes (r B /e) 2 
times longer. Using eq. (2), the time-scale becomes 4N 2/I3 
times larger than proposed by Moore et al. (1996). 

But there is another change to be made, as it is in 
fact the close encounters of halo particles that do the most 
damage to the subhaloes. According to Binney & Tremaine 
(1987), a good estimate for the disruption timescale can be 
had from an interpolation between the approximations for 
tidal encounters and for penetrating (b = 0) encounters. For 
each tidal encounter (Binney & Tremaine 1987, their eq. 
7-55), 
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while for each penetrating encounter (Binney & Tremaine 
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1987, their eq. 7-57) 



(AE) pcn = 



47rG m p ' 



;T, s (R)dR 



(11) 



% 7 (« 2 + ^ 

If the perturbed subhalo is an isothermal sphere, i.e 
E S (R) = v 2 /(6GR) « 0.2m h /(r hJ R), we find 



(AE) pc 



2G 2 m p m s r^ 1 



(12) 



Interpolating contributions from the tidal and penetrating 
encounters, i.e. 



AE= 4G 2 m 2 p m a r 2 



6 4 + ferf 



(13) 



finally allows us to integrate over all encounters. Following 
the procedure of Binney & Tremaine (1987), we simply find 
eq. (8) with r p (= e) replaced by 0.52(er^) 1/2 . Following 
Moore et al. (1996) again we get the same functional form 
as eq. (9), but the timescale is approximately 2(r s /e) 3//2 ~ 
times longer than estimated by Moore et al. (1996). 
By definition, m s /m p = N s , and we use eq. (2) to get 



t dis w 4.5iV s 6 — — « 2.2A" a 2 ^-e' 1 



(14) 



r s v h v h 
As ?"h is at least a few times r a , the disruption time is at least 
20-/V a 5//6 embedding halo crossing times. It is also a factor of 
18A r a _1 ' /6 In JVh ~ 100 times longer than the particle-subhalo 
two-body disruption timescale. 



4 CONCLUSIONS AND DISCUSSION 

Overmerging is the numerical disruption of subhaloes within 
embedding haloes. Of the three main two-body disruption 
processes, particle-subhalo two-body heating is clearly iden- 
tified as the cause for overmerging. Its timescale is shown to 
be much shorter than that for the two other processes, two- 
body evaporation and particle-subhalo tidal heating. Note 
that softened particles form into more extended subhaloes 
than is realistic, so they are more vulnerable to these disrup- 
tion processes (van Kampen 1995), and to possible physical 
disruption processes as well. 

Recently several research groups used simulations with 
a very high resolution in order to resolve the overmerging 
problem (Klypin et al. 1999a; Ghinga et al. 1998, 1999; 
Moore et al. 1999). Unfortunately, different group finders 
and different definitions for disruption times were used, so 
a direct comparison of the results is not straightforward. 
Still, the consensus is that increasing the number of parti- 
cles overcomes, at least partially, the overmerging problem. 
However, the resolution needs to be rather high: for N-body 
simulations on a cosmological scale, this requires the use 
of at least 10 9 particles, which is not very practical. Fur- 
thermore, for the smallest groups the overmerging problem 
simply remains. 

Another option is to include a baryonic component. 
With the addition of dissipative particles, haloes should be 
more compact and have a higher central density for the same 
numbre of particles. However, as Klypin et al. (1999a) re- 
mark, there is a limit to this as some fraction of the baryons 
tend to end up in rotationally supported disks. A more prac- 
tical problem with dissipative particles is the actual sim- 
ulation techniques needed, which usually is some form of 



smoothed particle hydrodynamics (SPH). The resolution of 
SPH codes is typically not as high as that of N-body codes, 
so for the purpose of resolving the overmerging problem it 
is not a useful alternative at present. 

A third option is to use halo particles, which prevents 
overmerging by construction (van Kampen 1995). The idea 
is that a group of particles that has collapsed into a viri- 
alised system is replaced by a single halo particle. Local 
density percolation, also called adaptive friends-of-friends, 
is adopted for finding the groups. This is designed to iden- 
tify the embedded haloes that the traditional percolation 
group finder links up with their parent halo. By applying 
the algorithm several times during the evolution, merging of 
already- formed galaxy haloes is taken into account as well. 
Once a halo particle is formed, more N-body particles will 
group around it at later times. If such a group can virialize, 
it is replaced by a more massive halo particle. This will usu- 
ally happen in the field. However, for halo particles that end 
up in overdense regions, the particles that swarm around a 
halo particles will be stripped quite rapidly. 

Once the overmerging problem is resolved down to the 
subhalo mass-scale one is interested in, the physical pro- 
cesses can be studied. This is becoming feasible for current 
simulations. However, whether the physical processes them- 
selves are properly modelled using N-body simulations has 
yet to be proven. The problem of artificially large subhaloes 
due to softening needs to be solved, for example. Another 
problem might be the modelling of dynamical friction, which 
requires a very smooth distribution of particles in the em- 
bedding halo in order to produce the wake that generates 
the drag force. 
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